Properties of the superconducting state in molecular metallic hydrogen under 

pressure at 347 GPa 
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The thermodynamic properties of the superconducting state induced in metallic molecular hydro- 
gen under the influence of pressure 347 GPa were determined. In particular, it has been shown that 
the critical temperature {Tc) changes in the range from 120 K to 90 K for fi* G (0.08, 0.15), where 
fj.* is the value of the Coulomb pseudopotential. Next, the energy gap near the temperature of 
zero Kelvin (2A (0)) was calculated. It has been stated, that the dimensionless ratio 2A (0) /fcsTcr 
slightly decreases with the increase of fj,* from 3.98 to 3.84. In the last step, the ratio of effective 
electron mass {ml) to the bare electron mass (me)) was determined. It has been proved that ml/me 
takes its highest value equal to 1.96 for T = Tc- 

PACS numbers: 74.20.Fg, 74.25. Bt, 74.62.Fj 



I. INTRODUCTION 

The influence of the pressure {p) on thermodynamic 
properties of the superconducting state is one of more 
interesting issues in the solid state physics. 

In the case of simple-metal superconductors with the 
electron-phonon pairing mechanism for a long time there 
was conviction, that the increase of the pressure leads 
to the decrease of the critical temperature (Tc) and 
the remaining thermodynamicparameters do not demon- 
strate interesting properties [l|. The change took place 
in the last years when in some of the classical materials 
the pressure-induced superconducting state of relatively 
high critical temperature was discovered, for example: in 
lithium [Tclfuax — 14 K for p ~ 30 GPa and in calcium, 
where [Tc]^^^ ~ 25 K for p ~ 160 GPa [2], [3]. 

The obtained results are keeping an eye on the metal- 
lic hydrogen again, where the existence of the pressure- 
induced superconducting state with the very high value 
of the critical temperature (comparable with the value 
of the room temperature) was expected [j|. From the 
theoretical point of view the following facts enable to 
make the above assumption true [5] : (i) high value of the 
phonon frequency caused by a small mass of atomic nu- 
clei building the crystal lattice (single proton), (ii) large 
value of the electron-phonon coupling related to the lack 
of the inner electronic shells and (iii) low value of the 
Coulomb pseudopotentiaHat least in the range of pres- 
sures up to 500 GPa 0, H)- 

Since the time of publishing the Ashcroft's work it is 
stated that the value of Tq for metallic hydrogen is lo- 
cated in the range from 80 K to 300 K for selected pres- 
sures lesser than 500 GPa [6|, 0, [3]. Under the ex- 
tremely high pressure {p = 2000 GPa) the metallic state 
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of the hydrogen may transit itself into the superconduct- 
ing state in the range of the temperatures from 413 K to 
631 K 0, [3]. Let us note that in this case the thermo- 
dynamic properties of the superconducting state strongly 
deviate from the predictions of the BCS theory [lO]. 

In the range of the "low" pressures (up to 500 GPa) 
the thermodynamic properties outside the critical tem- 
perature were not studied. For that reason we will cal- 
culate the selected thermodynamic parameters of the su- 
perconducting state in the molecular metallic hydrogen 
for the exemplary pressure of 347 GPa. In particular, 
in the framework of the Eliashberg approach [ll|, we 
will determine the value of the critical temperature as a 
function of the Coulomb pseudopotential. Next, we will 
exactly calculate the order parameter (A (T)) and the 
wave function renormalization factor. On the basis of 
obtained results the ratio 2A (0) /fesTc as a function of 
the Coulomb pseudopotential and the maximum value of 
the electron effective mass will be determined. Addition- 
ally, in the paper, the relationship between the structure 
of the electron-phonon coupling and the form of the solu- 
tions to the Eliashberg equations along the real axis will 
be discussed. 



II. THE ELIASHBERG EQUATIONS 

The Eliashberg equations can be determined on the 
real axis, imaginary axis or in the mixed representation 
(simultaneously on the real and imaginary axis) . Each of 
the mentioned approaches have the unique advantages, 
but also come with the characteristic mathematical prob- 
lems. The exact solution of the Eliashberg equations on 
the real axis composes an extremely difficult mathemat- 
ical problem [l2|. From the other side, the Eliahberg 
equations on the imaginary axis can be solved in a sim- 
pler way [l3|. However, in this case in order to interpret 
physically the obtained results (e.g. to calculate the value 
of the energy gap near the temperature of zero Kelvin) 



the achieved solutions should be analytically continued 
on the real axis. Unfortunately, the procedure of the an- 
alytical continuation is complicated and it demands very 
high precision during the numerical calculations '14']. It 
is worth mentioning, that the results obtained in a such 
way are stable only in the range of the low frequencies. 

The Eliashberg equations in the mixed representation 
are proclaimed to be a reasonable compromise between 
the two previously mentioned approaches. Firstly, they 
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can be exactly solved in the much simpler way than equa- 
tions on the real axis. Secondly, the obtained results are 
stable even for the very large values of the frequency. In 
the case of molecular metallic hydrogen it is an important 
matter, because the maximum phonon frequency (r^max) 
is equal to 477 meV 0. 

The Eliashberg equations in the mixed representation 
can be written in the form 11 ll , llSlI : 
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In Eqs. (HI and ^ the symbols (j) and Z denote the 
order parameter function and the wave function renor- 
malization factor determined respectively on the real axis 
{uj) or imaginary axis {uj„i), where aj,„ = (7r//3) {2m — 1) 
is the TO-th Matsubara frequency and /3 = {ksT) 
{ks is the Boltzmann constant). The order parameter 
on the real axis (A {uj)) is defined in a following way: 
A{uj) = (f) {uj) / Z {oj) . The pairing kernel of the electron- 
phonon interaction is given by the expression: 
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where the Eliashberg function {a^F (51)) for the molec- 
ular metallic hydrogen at the pressure of 347 GPa was 
determined in the paper [9| . Symbol /i* denotes the value 
of the Coulomb pseudopotential; 9 is the Heaviside unit 
function and uJc is called the phonon cut-off frequency 
{ujc — 351,iiax). In the presented work, /z* was taken into 
account parametrically. In particular, it was assumed 



that: ^* e (0.08,0.15). In the last two terms of Eqs. (H]) 
and ([5]) the symbols N {oj) and / {uj) denote the statisti- 
cal functions of bosons and fermions respectively. 

The Eliashberg equations were solved numerically for 
M = 800. In this case the stability of all solutions was 
assured for the minimal temperature equal to 23.2 K (2 
meV) regardless of the assumed value of /i* . 



III. THE NUMERICAL AND ANALYTICAL 
RESULTS 

A. Value of the critical temperature 

The dependence of the critical temperature on the 
Coulomb pseudopotential for the low-Tc superconduc- 
tors can be determined in the easiest way when using the 
formula drawn out by Allen and Dynes [16|. Unfortu- 
nately, in the case of the metallic hydrogen, the classical 
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FIG. 1: The value of the critical temperature as a function 
of the Coulomb pseudopotential. The filled circles repre- 
sent the exact results obtained with an use of the Eliashberg 
equations. The solid line has been obtained using the modi- 
fied Allen-Dynes formula. Empty squares and triangles refer 
respectively to the classical Allen-Dynes expression and the 
McMillan formula (/i = /2 = 1) [l^, [l7|. 



formula determines correctly Tc {^J-*) only for the very 
low values of fi* (see Fig.[l]); for the highest values of 
the pseudopotential the critical temperature is consider- 
ably underestimated. For that reason selected parame- 
ters that appear in a classical expression were calculated 
once again. In particular, we have used the least-squares 
analysis and 120 values of Tc (/i*) obtained with a help 
of the Eliashberg equations. In the considered case, the 
modified Allen-Dynes formula takes the form: 
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where the strong-coupling correction function (/i) and 
the shape correction function (/2) are given by: 
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The new functions Ai and A2 can be written as follows: 
Ai = 4.11(l-4.6lAi*) (7) 

and 
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The parameter W2 denotes the second moment of the nor- 
malized weight function and should be calculated on the 



basis of the formula: 
2 
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The quantity uji^ is called the logarithmic phonon fre- 
quency and A is the electron-phonon coupling constant: 
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In the case of molecular metallic hydrogen {p — 347 GPa) 
the following results were achieved: -^0^2 = 199 mcV, 
Win = 142 meV and A = 0.927. 

Below, we draw the reader's attention to the fact, that 
the new parameterization of Ai and A2 very relevantly 
affects the dependence of the functions /i and /2 on the 
value of the Coulomb pseudopotential. Namely, it comes 
to: fi > f2 — 1. Additionally, the values of /i grow 
together with the growth of the pseudopotential. The 
above result means that Tc very significantly depends 
on the strong-coupling effects modeled by the function 
/i ; while the shape correction function weakly affects the 
critical temperature. 

When coming back to the data presented in Fig. [1] it 
can be easily noticed, that in the considered range of the 
Coulomb pseudopotential' s values the critical tempera- 
ture decreases from 120 K to 90 K. The achieved result 
clearly states that, even for the relatively large /i*, the 
critical temperature takes very high value. 



B. The order parameter function 

The Eliashberg equations were solved for the range of 
temperature from 23.2 K to Tc- As an example, the 
dependence of the real and imaginary part of the order 
parameter on the frequency in Fig. [5] was presented. For 
the low values of oj the non-zero is only the real part of 
the order parameter. Next, we observe the characteris- 
tic sequence of the local maximums (minimums) and the 
points of the fold both for Re[A (w)] and Im[A (w)]. Ad- 
ditionally, in Fig. [5] the Eliashberg function is shown; in 
order of the better representation its values were multi- 
plied by 15. When comparing the shapes of Re [A (w)] 
and Im[A (w)] with the form of the Eliashberg function 
the correlation between the characteristic points of the 
considered functions can be easily noticed. Let us mark 
the fact, that from the physical point of view, the exis- 
tence of the peak in the Eliashberg function determines 
the area of to in which the electron-phonon interaction is 
exceptionally strong. Thus, the distribution of the char- 
acteristic points of the real and imaginary part of the 
order parameter is inseparably connected with the spe- 
cific structure of the electron-phonon interaction. 



35 

30 

25 

20 

15 

10 

5 



-5 

-10 

-15 




■Re[A(co)] meV 
■ lm[A(a))] meV 

15a'F(co) 







100 



200 300 

CO meV 



400 



-TU 


□ T=23.2K 


35 


/^"^ o T=81.2K 
/ \ ^ T=98.6K 


30 


/ ^ O T=104.4K 
" / ? ft T^=108.2K 


^25 


- /-"xl 


^20 




/ 


A^ 


i15 
E 10 




I 




5 


w 


hii 





- \^.^a 



-30 -15 15 30 
Re[A((o)] meV 



18 


t"-—,,^ 






16 


: 




• 


14 


• 
• 




• 


>12 


• 


CD 


• 


E10 


• 


^— V 


• 


b 8 


• 


< 


• 


6 


• 


4 


• 


2 

n 


, . M: 



20 40 60 80 
TK 



100 



FIG. 2: The real and imaginary part of the order parameter 
as a function of the frequency for T = 58 K (5 meV) and 
fj,* = 0.1. In the figure there is also plotted the rescaled 
Eliashberg function. 



In a global way the dependence of A (w) on the temper- 
ature was presented in Fig. [3] (A). It can be clearly seen 
that on the complex plane the values of A (w) form the 
characteristic "ear". The obtained result allows to char- 
acterize the effective potential of the electron-electron in- 
teraction, which is connected with the real part of the or- 
der parameter [18|. In particular, let us notice that only 
in the range of frequencies from zero to the frequency 
slightly lesser than flmax, the effective electron-electron 
interaction is attractive (Re[A (w)] > 0); the achieved 
conclusion is true for any considered value of fi* . 

Below, the values of the ratio i?i = 2A (0) /ksTc for 
fi* e (0.08,0.15) are determined; the symbol A (0) de- 
notes the order parameter near the temperature of zero 
Kelvin. Let us notice, that in the framework of the weak- 
coupling approach (the BCS model) the parameter Ri 
takes the constant value equal to 3.53. In the case when 
A > 0.2 its exact value can be calculated only with the 
help of the Eliashberg equations. In particular, A (T) is 
calculated on the basis of Re [A (w)] with an use of the 
expression; 



A (T) = Re [A (tj = A (T))] 



(12) 



In Fig. [3] (B) the dependence of the order parameter on 
T is plotted. Due to the fast saturation of the function 
A (T) in the area of the lower temperatures, for A (0) 
one can assume the value of A (T = 23. 2K). The final 
results are presented in Fig.|4l where the function i?i (fi*) 
is shown. It is easy to notice, that the values of the ratio 
Ri are higher than [i^ilg^g and only slightly decreasing 
with the growth of the Coulomb pseudopotential (from 
3.98 to 3.84). The determined dependency of i?i on fi* 
is connected with the strong, but comparable, influence 
of the electronic depairing correlations on Tc and A (0); 
see Fig. [T] and the inset in Fig. SI 



FIG. 3: (A) The order parameter on the complex plane for the 
selected values of the temperature. The solid lines with points 
were obtained for cj e (0, fimax); the solid lines represent the 
results for the frequency from f^max to 1600 meV. (B) The 
dependence of the order parameter on the temperature. In 
both cases /i* =0.1 was assumed. 
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FIG. 4: Dimensionless ratio Ri as a function of /i*. The de- 
pendence of the parameter A (0) on the Coulomb pseudopo- 
tential has been shown in the figure's inset. 



C. The electron effective mass 

In Fig. [5] the wave function renormalization factor is 
plotted. Similarly as for the order parameter, in the low- 
frequencies region the non-zero is only the real part of 
Z (oj). For higher frequencies, we can observe the com- 
plicated dependence of Re[Z] and Im[2'] on oj, which is 
plainly correlated with the form of the Eliashberg func- 
tion. However, it should be noticed, that in the opposi- 
tion to the order parameter, the function Z (cu) is signif- 
icantly weaker dependent on the temperature (see Fig. 

In the framework of the Eliashberg theory the wave 




200 300 

CO meV 



FIG. 5: The real and imaginary part of the wave function 
renormalization factor as a function of frequency for T = 58 
K and fj* = 0.1. In the figure the Eliashberg function is also 
plotted. 
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FIG. 6: The wave function renormalization factor on the com- 
plex plane for the selected values of T. In the inset the de- 
pendence of ml /me on the temperature has been marked. In 
both cases fi* = 0.1 was assumed. 



fective mass in the superconducting area the dependence 
of ml/m,, on the temperature is plotted in the inset of 
Fig. [6] For the molecular metallic hydrogen the value of 
the electron effective mass is pretty high and takes its 
nraximum equal to 1.96 for T = Tc- Let us notice that 
in this case ttt.* is independent of /i* . 
IV. SUMMARY 



In this paper the thermodynamic properties of the su- 
perconducting state induced in molecular metallic hy- 
drogen (p = 347 GPa) were studied. In order to do 
that, the Eliashberg equations in the mixed represen- 
tation were exactly solved. The obtained results en- 
abled to determine the dependence of the order param- 
eter and the wave function renormalization factor in the 
wide range of the frequency. It has been stated, that for 
w e (0, f^niax) the existence of the characteristic points of 
the function A (w) and Z (uj) is strictly correlated with 
the structure of the electron-phonon coupling modeled 
by the Eliashberg function. Additionally, in the range 
of frequencies from zero to the frequency slightly lesser 
than llmax, the effective electron-electron interaction is 
attractive (Re[A (w)] > 0). 

When basing on the exact solutions of the Eliashberg 
equations it has been shown, that the value of the criti- 
cal temperature decreases slower with the growth of the 
Coulomb pseudopotential, than it is predicted by the 
classical Allen-Dynes formula. In particular, Tc changes 
from 120 K to 90 K for n* G (0.08,015). Next, the de- 
pendence of the dimensionless ratio Ri on /u* was deter- 
mined. It has been proven that the parameter i?i weakly 
depends on the value of the Coulomb pseudopotential 
and [Ri],^^^ = 3.98 for fi* = 0.08. In the last step, the 
value of the electron effective mass in the superconduct- 
ing area was calculated. On the base of the obtained 
results it has been shown, that [m-e/"^e]max ~ ^-^^ ^'^^' 
T = Tc. 
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function renormalization factor plays a very significant 
role, because with the help of it one can calculate the ra- 
tio of the electron effective mass (m.*) to the bare electron 
mass (we): 
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Re [Z (0)] . 



(13) 



In order to precisely track the values of the electron ef- 
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